Mon. Not. R. Astron. Soc. 000,11]-?? (2006) Printed 5 February 2008 (MN WT^ style file v2.2) 



Spinning-Down of Moving Magnetars in the Propeller 
Regime 



o 
o 

(N 
=3 



O.D. Toropina/*, M.M. Romanova,^, and R.V.E. Lovelace, ^'^ 

^ Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, Moscow 117997, Russia 
^ Department of Astronomy, Cornell University, Ithaca, NY 14853-6801 

^ Department of Applied and Engineering Physics, Cornell University, Rhaca, NY 14853-6801 



5 February 2008 



(N 



Key words: neutron stars — magnetars — magnetic field — MHD 



> 

(N 

o 

\o 
o 

o 



ABSTRACT 

We use axisymmetric magnetohydrodynamic simulations to investigate the 
spinning-down of magnetars rotating in the propeller regime and moving supersoni- 
cally through the interstellar medium. The simulations indicate that magnetars spin- 
down rapidly due to this interaction, faster than for the case of a non-moving star. 
From many simulation runs we have derived an approximate scaling laws for the an- 
gular momentum loss rate, L oc — 7]^^ ffi-^ p^-^ j\4~'^ '^nl-^ , where p is the density of the 
interstellar medium, A4 is Mach number, p is the star's magnetic moment, 51* is its 
angular velocity, and rjm is magnetic diffusivity. A magnetar with a surface magnetic 
field of 10" - 10" G is found to spin-down to a period P > lO'^ - 10*^ s in ~ Iff* - 10^ 
years. There is however uncertainty about the value of the magnetic diffusivity so 
that the time-scale may be longer. We discuss this model in respect of Soft Gamma 
Repeaters (SGRs) and the isolated neutron star candidate RXJ1856. 5-3754. 



1 INTRODUCTION 



Some neutron stars referred to as "magnetars" have unusu- 
ally large magnetic fields, B ~ lO" - 10^^ G (Duncan & 
Thompson 1992; Thompson & Duncan 1995). Possible can- 
didates for magnetars include anomalous X-ray pulsars and 
soft gamma-ray repeaters (SGRs) (Kulkarni & Frail 1993; 
Kouveliotou et al. 1994; Hurley et al. 1999). These objects 
are associated with supernovae remnants and hence are rel- 
atively young (Vasist & Gotthelf 1997; Kouveliotou et al. 
1998). Only a few candidates for magnetars have been found 
so far. The estimated birthrate of magnetars is ^ 10% of or- 
dinary pulsars (Kulkarni & Frail 1993; Kouveliotou et al. 
1994, 1999) so that there might be many more magnetars 
which are presently invisible. Their "visibility" depends on 
a number of factors. One important factor is the rate of the 
star's spin-down. If magnetars spin-down rapidly to very 
long periods, then one will not detect spin modulated vari- 
ability during fiares. Consequently the identification of the 
variability with a rotating neutron star would be more dif- 
ficult. 

During the pulsar stage of evolution, magnetars spin 
down much more rapidly than ordinary pulsars. Conse- 
quently, they pass through their pulsar stage much faster, 
in ~ 10* years (Thompson & Duncan 1995). When the 
light cylinder radius becomes larger than magnetospheric 



radius Vm, the relativistic wind is suppressed by the in- 
flowing matter (Shvartsman 1970) and the star enters the 
propeller regime where the spin-down is due to the interac- 
tion of the star's rotating field with the interstellar medium 
(ISM) (Davidson & Ostriker 1973; lUarionov & Sunyaev 
1975). Magnetars with velocities w > 100 - 200 km/s in- 
teract directly with the ISM. That is, the magnetospheric 
radius is larger than gravitational capture radius (Harding 
& Leventhal 1992; Rutledge 2001; Toropina et al. 2001). In 
the propeller regime the rapidly rotating magnetosphere in- 
teracts strongly with the supersonically inflowing ISM. 

The spin-down rate of supersonically moving magne- 
tars in the propeller regime has been estimated earlier, but 
different authors have obtained rather different results. For 
example, Rutledge (2001) estimates a spin-down time of 
~ 4 X 10^ yr for a neutron star with a surface magnetic 
fleld B = lO^'^ G and velocity v = 300 km/s. On the other 
hand, Mori and Ruderman (2003) estimate that a magnetar 
spins-down to periods greater than 10* s within ~ 5 x 10^ 
years. Mori and Ruderman put forward this model as an ex- 
planation of the isolated neutron star (INS) candidate RX 
J1856. 5-3754 which does not show variability. 

The propeller stage of evolution has been investigated 
in non- magnetar cases both theoretically (e.g., lUarionov 
& Sunyaev 1975; Davies, Fabian & Pringle 1979; Davies 
& Pringle 1981; Lovelace, Romanova & Bisnovatyi-Kogan 



2 O.D. Toropina et al. 





, , , , I , . . . . 1^ 1^ . . . J . 1^ . 1^ I I I k^^^H 

-0.5 0.5 2 1 1-5 ^°gioP 

Figure 1. Matter flow around a strongly magnetized star rotating in the propeller regime and propagating through the interstellar 
medium with Mach numbers Ai = 1 (top panel) and M = 3 (bottom panel). Other parameters correspond to the main case. The 
background represents the logarithm of density and the length of the arrows is proportional to the poloidal velocity. The solid lines 
are magnetic field lines. The dashed solid line shows the Alfven surface. Distances are measured in units of the Bondi radius. Time 
corresponds to 50 rotation periods of of the star. 



1999; Ikhsanov 2002; Rappaport, Fregeau, & Spruit, 2004) 
and with MHD simulations (Wang & Robertson 1985; Ro- 
manova et al. 2003; Romanova et al. 2004; Romanova et al. 
2005; Ustyugova et al. 2006). However, only limited theoret- 
ical work has been done for magnetar-type propellers, which 
propagate through the ISM supersonically (Rutledge 2001; 
Mori & Ruderman 2003). No MHD simulations of this flow 
regime have been done previously. 

This work presents the first numerical simulations of 
the interaction of supersonic, fast rotating magnetars with 
the ISM. Earlier, we investigated supersonic propagation 
of non-rotating strongly magnetized stars through the ISM 
(Toropina et al. 2001). The present simulations are analo- 
gous to those in Toropina et al. (2001). However, here the 
star rotates in the propeller regime. The main objective of 
this work is to determine the dependences of spin-down rate 
of the star on the different variables and to estimate corre- 
sponding time-scale of the spin-down. 

Sections 2 and 3 describe the physical situation and 



simulation model. In §4, we discuss the results of our simu- 
lations. In §5, we apply our results to magnetars and in §6 
we discuss possible magnetar candidates. Conclusions are 
given in Section 7. 



2 PHYSICS OF THE "PROPELLER" REGIME 
OF A MOVING MAGNETAR 

Here, we consider the spin-down of a highly magnetized ro- 
tating neutron star or "magnetar" moving at a velocity v 
faster than the local sound speed of the interstellar medium; 
that is, 11 > Cs, where Cs ~ 10*' cm/s for a 10* K hydrogen 
plasma. Theory and simulations of the spin-down for the 
case of Bondi accretion to a non-moving, rotating magne- 
tized star in the propeller regime was discussed earlier by 
Romanova et al. (2003). A general discussion and review 
of theory of the spin-down of moving rotating magnetars is 
given by Mori and Ruderman (2003). 
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Figure 2. Temporal variation of the total angular momentum flux from the star obtained by integrating over a surface surrounding the 
star (solid line) and the flux through a surface a<;ross the magnetotail at z = 0.6 (dotted line). 



We first discuss tlie diflFerent lengtii scales in tiiis prob- 
lem. We let Tm denote the ciiaracteristic size of tlie star's 
magnetosphere witli the region inside rm magnetically dom- 
inated and corotating with the star at angular rate f2* . 

The standoff distance of the bow shock formed by the 
interaction of the moving ISM with the star's magnetic field 
is at a distance of the order of 

from the center of the star. Here, /X33 = p/lO'^'^Gcm^ is the 
star's magnetic moment, n is the number density of the ISM 
in H/cm^, and vy = v/(10'^cm/s). Note that the magnetic 
field at the star's surface is B = fi/R^ = 10^V33(10km/i?.)^ 
G, where i?* is the star's radius. We assume that Vst < r-m- 
In the propeller regime the magnetospheric radius Vm 
is larger than the corotation radius 

r„ = {GM/^tf^ ^ 1.7 X 10"' Pg''^ cm , (2) 

where P is the rotation period of the star and P3 = 
s. We assume M = IAMq here and subsequently. Conse- 
quently, the incoming matter is flung away from the star in 
the equatorial plane. 

A further characteristic length is 

= v/n, = 1.5 X 10^7^:^ cm , (3) 

This is the distance the star moves in a rotation period 
divided by 2n. A related length is 

Vs = Cs/0,* , (4) 

where Cs is the sound speed in ambient ISM. Other char- 
acteristic lengths can be constructed from combinations of 
equations (1) - (4). For example, the Bondi radius is given by 
re = fcr/^'s, and the Bondi-Hoyle radius by rsH = r^r/r^- 

Additionally, there is a length scale associated with the 
magnetic diffusivity 77 

rr, = . (5) 



Roughly, is the distance the magnetic field can diffuse 
in a rotation period of the star. These lengths are expected 
to be small compared with the lengths (1) - (4). Note that 
magnetic diffusivity is much more important in this problem, 
compared to viscosity, so that we do not take into account 
viscosity in the analysis below. 

Note also that wc assume the light cylinder radius, 

rL = cP/27r ^ 4.8 x lO^^Pa cm , (6) 

is larger than the standoff distance Vst- 



2.1 Torque Derived from Dimensional Analysis 

The rate of spin-down of the star is estimated as dL/dt = 
IdQt/dt Si — 47rr^/rmBm/47r, which is the magnetic torque 
at the radius rm, where / is the star's moment of inertia, 
and / is the fraction of the solid angle where this torque 
acts. We have assumed that « Bm at rm- Thus 

dL^^dn^^_fjl_ 

dt dt ^ ' ^ ' 

We take an empirical approach to determining the depen- 
dences of the magnetospheric radius Vm on the different 
physical parameters 

Tm =rm{p, M. Cs, O, , Tj) , (8) 

as suggested by Mori and Ruderman (2003). We write 

rm = {r.tr{r,rf{rr,y{rsr{r,y , (9) 
where s = l — a — f3 — f — e. Consequently, 

^ oc - p"/V'""r?"'^^'c7^^t;^=-"a2 , (10) 

where g = 2/3 -|- 37/2 -|- 3e -|- 3s. In §4 we discuss the values 
of the exponents suggested by our simulations. 
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Figure 3. Distribution of the angular momentum fluxes in thie magnetotail for our main case. The color background shows the specific 
angular momentum carried by the magnetic field (top panel) and by the matter (bottom panel). The solid lines are magnetic field lines. 



2.2 Torque Derived from Analytical Models 

There is a wide variety of theories which predict different de- 
pendencies for the spin-down torque. The dependencies are 
summarized in tables by Davies, et al. (1979) and by MR03. 
For example, from the table in MR03 we see that the power 
law dependence on varies strongly, from no = —1 up to 2; 
the power law dependence on ^ varies less dramatically from 
=2/3 up to 5/3. The dependence on the star's velocity 
varies dramatically from n„ = —5/3 to 7/3. The dependence 
on the density of the ISM varies from np = 1/6 to 2/3. The 
torque obtained from the different models differs by several 
orders of magnitude. The spin-down time scales vary cor- 
respondingly. One of the reasons for the differences is that 
in different situations (disk accretion, wind accretion to the 
non-moving or moving pulsar) the physics of interaction is 
different. Only a few papers deal with the situation consid- 
ered in this paper where the magnetospheric radius is larger 
than the gravitational capture radius. One consequences of 
this is that the magnetospheric radius Tm is of the order of 
the stand-off radius rst ■ Below, we propose a simple physical 
theory which is closest to our particular problem. 



2.3 Torque Derived from Simple Physical Model 

Ifere, we give a simple model for the spin down of a moving 
propeller star. We assume that the magnetospheric radius is 
larger than the accretion radius so that the magnetosphere 
is an obstacle for the interstellar matter. Matter stops at 
the stand-off distance rat determined from the balance of 
the ram pressure of the ISM and the magnetic pressure of 
the star and is given by equation (1). The mass flow crossing 
the area vrr?, is 



M — -Kr^tpv 



We assume that a fraction / = k[ri/{rstv)Y of this mass 
flux participate in spinning-down of the star, where Rem = 
rstv/rj is a magnetic Reynolds number and k and s are con- 
stants. The factor / is determined by the penetration of 
the interstellar plasma into the star's magnetosphere. The 
larger the penetration the larger the transfer of angular mo- 
mentum from the star to the ISM. Noting that Q*rat is the 
specific angular momentum of the matter at a distance rst 
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Figure 4. Dependence of the angular momentum flux on difl'erent parameters, (a) the Mach number M, (b) the angular velocity of the 
star ujt = f2*/r2x, (c) the magnetic moment fi, and (d) the magnetic diffusivity ijm- 



we obtain 



3 NUMERICAL MODEL 



— = -fc — A'm,rlt~~krj''p—^a^-v 5 n^, 

at \rstvJ 

As an illustration consider s = 0.3 which is suggested by our 
simulations (see Figure 4d). We then find 



(11) 



This model indicates that the torque increases with /i, p and 
fl*, but it but decreases with increasing velocity v. 



2.4 Transition from Pulsar to Propeller Stages 

Equating the standoff distance Tst to the fight cylinder ra- 
dius rL gives an estimate of the transition from the pulsar 
stage to the propeller stage, 



Ptrans = 5.6 X 10 



1/3 
2 / M33 



1/fi 1/3 



(12) 



The propeller stage will end when the corotation radius Vcr 
is larger than the magnetospheric radius Tst ■ This occurs for 



P = 2.1 X 10 



1/2 

,6 / M33 



1 1/2 



(13) 



Thus, a magnetar with field B — 10^^ G propagating 
through the ISM supersonically will be in the propeller 
regime for rotation periods 560s < P < 2 x lO^s. Our model 
applies only for this range. 



We investigate the interaction of fast moving rotating mag- 
netized star with the ISM using an axisymmetric, resistive 
MHD code. The code incorporates the methods of local it- 
erations and fiux-corrected-transport (Zhukov, Zabrodin, & 
Feodoritova 1993). The flow is described by the resistive 
MHD equations: 



dp 
dt 



-KV-(p v) = , 



P^+P(V- V)V: 



dt 



-Vp+ - J X B + F'' 
c 



f ^Vx(vxB) + ^V^B, 

^ + V.(psv) = -pV.v+^! 
at a 



(14) 



We assume axisymmetry {d/d(j} — 0), but calculate all 
three components of velocity and magnetic field v and B. 
The equation of state is taken to be that for an ideal gas, 
p = (7 — l)p£, with specific heat ratio 7 = 5/3. The equa- 
tions incorporate Ohm's law J = cr(E -I- v x B/c), where a 
is the electrical conductivity. The associated magnetic diffu- 
sivity, T]m = c^/(47r(j), is considered to be a constant within 
the computational region. In equation (2) the gravitational 
force, = —GMpH/R^ , is due to the central star. 

We use a cylindrical, inertial coordinate system (r, tj), z) 
with the 2— axis parallel to the star's dipole moment fl and 
rotation axis fl. The vector potential A is calculated so that 
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V-B = at all times. The star rotates with angular velocity 
f2» = f2» z. The intrinsic magnetic field of the star is taken 
to be an aligned dipole, with vector potential A = /xxR/ii^. 

We measure length in units of the Bondi radius Rb = 
GAf/Cj, with Cs the sound speed at infinity, density in units 
of the density of the interstellar medium p, and magnetic 
field strength in units of Bo which is the field at the pole 
of the numerical star. The magnetic moment is measured in 
units of ^0 = B(]R%/2. 

After reduction to dimcnsioiilcss form, the MHD equa- 
tions involve dimcnsionless parameters: 



SttPq 



RbVo RCr, 



(15) 



where fim is the dimcnsionless magnetic difi^usivity, the Rem 
is the magnetic Reynolds number, and /3 is so-called gravi- 
magnetic parameter, which describes the ratio between pres- 
sure of the ISM medium and magnetic field at the poles of 
the "numerical" star, B[). 

Simulations were done in a cylindrical region {Zmin ^ 
z < Zmax,0 < r ^ Rmax)- The numerical star was repre- 
sented by a small cylindrical box with dimensions R* « 
Rmax and \Zt \ « Zmax ■ A Uniform (r, z) grid with Nr x Nz 
cells was used. 

Initially, the density p{r,z) and the velocity v(r, are 
taken to be constant in the region p{r,z) = p and w = «*, 
Vtj, = 0. Also, initially the vector potential A was taken to 
that of a dipole so that — 0. The vector potential was 
fixed inside the numerical star and at its surface during the 
simulations. The star was initialized to be rotating at the 
rate fi*. 

The outer boundaries of the computational region were 
treated as follows. Supersonic inflow with Mach number A4 
was specified at the upstream boundary {z = Zmin , < r ^ 
Rmax)- The variables (p, Vr, Vz, e) are fixed. The inflowing 
matter is unmagnetized with B = 0. At the downstream 
boundary {z — Zmax,0 ^ r ^ Rmax) and at the cylindrical 
boundary {Zmin ^ z ^ Zmax,r = Rmax a free boundary 
condition was applied, d/dn = 0. Boundary conditions are 
described in greater detail in TOl and R03. 

The size of the computational region for most of the 

simulations was Rmax = 0.9Rb — 0.9, Zmin = —Rmax = 

-0.9, with Zmax = 2Rmax = 1.8. The grid Nr x Nz was 
289 X 865 in most of CEises. The radius of the numerical star 
was Rt = 0.025i?B = 0.025 in all cases. 



4 RESULTS 

We investigated numerically interaction of the rotating mag- 
nctar with the ISM and the rate of spinning-down. Simula- 
tions were done at a variety of parameters: magnetic mo- 
ments of the star p, angular velocities f2* , Mach numbers 
Ai and diffusivities fj. In most cases we varied one param- 
eter in a time, and kept other parameters fixed and corre- 
sponding to the "main" case. In the main case a star ro- 
tates with an angular velocity ui, = i^t/^^K* ~ 0.7, where 
Qk* = \/ GM/Ri is Keplerian angular velocity at the sur- 
face of the numerical star. We suggest that numerical star 
(inner boundary) is much larger than the true radius of the 
neutron star. Mach number is = 3. We take gravimag- 
netic parameter f3 = 10~® which corresponds to the dimcn- 



sionless magnetic moment p, = 10 and take magnetic 
diffusivity fjm = 10~^. 



4.1 Dependence of Torque on Mach Number 

For parameters corresponding to our "main case," we did 
simulation runs for the Mach numbers M = 1, 3, 5, 6, and 
10, with the ambient sound speed fixed. From these runs 
we find that the torque decreases with the Mach number 
approximately as L oc —M~°''^. 

Figure la shows that for a Mach number M = 1 the 
flow is similar to that observed in case of a non-moving pro- 
peller (Romanova et al. 2003) . Namely, the rapidly rotating 
magnetosphcrc pushes matter and magnetic flux outward in 
the equatorial plane forming the low-density, rotating torus 
in the equatorial plane. The gravitational radius is several 
times larger than magnctosplicric radius so that a significant 
part of the inflowing matter is gravitationally trapped and 
accumulates around the star. This is similar to the case of 
spherical Bondi accretion to a star in the propeller stage. 
For M = 1, there is an axial flow of matter downstream 
from the shock wave. However, the cncrgy-deiisitj' of the in- 
flowing matter pv^/2 is smaller than the energy-density of 
the equatorial propeller outflow pv^/2 -|- B^/Stt and this is 
why the equatorial structure forms. 

Figure lb shows the flow at a larger Mach number 
= 3. In this case the energy-density of the ISM matter is 
laxger than the energy-density of the equatorial propeller- 
generated wind, so that the disk structure is bent and 
pushed to the direction of the tail. This interaction is similar 
to that observed in the simulations of magnetized supersonic 
stars in the non-rotating case (Toropina et al. 2001). Namely, 
the magnetosphere of the star acts as an obstacle for the ISM 
matter so that a bow shock stands in front of the star and a 
conical shock wave forms behind it. The stand-off distance 
and the cross-section of the interaction is larger in the pro- 
peller case than in the non-rotating case. This is because the 
rotating equatorial disk of matter and magnetic field gener- 
ated by the fast rotating magnetosphere. For Mach numbers 
M > 2 — 3, this disk structure is pushed by the inflowing 
matter into a magnetotail behind the star. The maximum 
rotational velocity of the magnetosphere Vprop = rD,* is sev- 
eral times larger than the velocity of the ISM matter, v. In 
common terminology the propeller is "supersonic" . In spite 
of this, the ISM matter pushes this fast rotating matter to- 
gether with the magnetic held into the magnetotail. Because 
of angular momentum conservation, the matter in the mag- 
netotail continues to rotate. In addition, the magnetic field 
is twisted by the rotating matter. The amount of the twist 
depends on the parameters. For the parameters of our main 
case the twist is B^/Bp ~ 0.1. 

For even larger Mach numbers, M = 6 — 10, the flow 
is similar to that observed for TVf = 3. However, the stand- 
off distance is even smaller, because the larger portion of 
the rotating magnetosphere is pushed into the magnetotail 
behind the star. 

The rotating magnetosphere interacts with the non- 
rotating matter of the ISM and this leads to the spinning- 
down of the star. Angular momentum lost by the star flows 
out from the star. Thus, we calculate the torque on the star 
by integrating the angular momentum flux density over a 
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surface surrounding the star, 

L = -|dS.(pv.™,-5^) . 

Here, dS is the outward pointing surface area element and 
the p— subscript indicates the poloidal component. Wo use 
a cylindrical surface around the star approximately at the 
Alfven surface to evaluate this integral. Figure 2 (solid line) 
shows temporal variation of this flux. One can see that 
after few initial rotations of the star, L becomes approxi- 
mately constant. For comparison, wc also calculated angu- 
lar momentum flux through the magnetotail at the distance 
z = 0.6 from the star. We obtained similar value of the 
angular momentum flux (see dashed line at the Figure 2) 
because the angular momentum lost by the star flows into 
the magnetotail. 

Figure 3 shows distribution of angular momentum flux 
densities carried by the magnetic field — BprB^/47r (top 
panel) and that carried by the matter pvprv^ (bottom 
panel). One can see that inside the star's magnetosphcrc 
the magnetic field gives the main contribution to the angu- 
lar momentum flux (see top panel). After passing the shock 
wave, magnetosphere interacts with non-rotating ISM and 
transports part of its angular momentum to the matter. This 
rotating matter propagating to the magnetotail gives the 
main contribution to the angular momentum flux in the tail 
(see bottom panel) . The rate of angular momentum loss from 
the star depends on efficiency of mixing of the ISM matter 
with the magnetic field of the magnetosphere. 



4.3 Summary of Scaling Laws 

Taking into account all above dependencies, we obtain the 
summary of scaling laws: 

^ oc - « . (16) 

The dcpondoncios do not exactly coincide with those 
derived in the simple physical model (§2.3, eq. 11), however 
the main tendencies do coincide. It is important to note that 
in both cases the dependence on Mach number (velocity of 
the star, v, in the physical model) is negative. This is a robust 
result which was understood analytically and confirmed in 
simulations. 

We also can compare these empirical dependencies with 
scaling laws proposed in §2.1. The dependence of L on the 
magnetic moment /x implies that a « 1.4. The dependence of 
L on the magnetic diff'usivity rj implies that 7 ~ —0.2. The 
dependence of L on velocity of the star v or Mach number 
A4 implies that /? -|- e « —0.8. The dependence of L on the 
star's rotation rate 0» implies that f3 « —2.4. Thus equation 
(10) implies that 

dL 0.7 0.6 0.3 -3e -0.4ol-5 ria\ 

— oc-purjc^v S2». (16) 
at 

Note, that dependence on density p"'^ was derived from the 
scaling analysis and it is very close to that obtained from 
the numerical simulations. 

The dependences we find of L on /j, v, and f2 are not 
compatible with the equation (4) of Mori and Ruderman 
(2003) because different approaches were used. 



4.2 Dependence of the Torque on Other 
Pcirameters 

We performed a number of simulations with different angu- 
lar velocities of the star fi, and different magnetic moments 
/^. We calculated the total angular momentum loss rate from 
the star and investigated its dependence on the 0, and /x. 

Figure 4b shows the total angular momentum loss rate 
from the star as a function of fi, . We find that for the main 
case (Mach number A4 = 3) this dependence is L oc —ftl'^. 
This dependence is somewhat stronger than that found for a 
non-moving star where L oc — fjj,'^ (Romanova et al. 2003). 
We expected that at Mach number M. = 1 the torque will 
have a power between 1.5 and 1.3, however, we obtained q = 
1.1 value (see Figure 4b, rhombs). This may be connected 
with different physics of interaction at very large and very 
small Mach numbers. 

Figure 4c shows the total angular momentum loss rate 
from the star as a function of p. Approximately, L oc —pP'^. 
This dependence is similar to that for the case of non-moving 
star L oc — /u"'^ (Romanova et al. 2003). 

We also calculated the dependence of the torque on dif- 
fusivity, and obtained: L oc — r;" Diffusivity is important 
factor which helps to mix non-rotating matter of the ISM 
with rotating magnetosphere. It is one of important fac- 
tors which determines the amount of matter participating 
in spinning-down. 

In addition we calculated the dependence of the torque 
on the density of the ISM medium and obtained dependence: 
L oc 



5 SPINNING DOWN OF MAGNETARS 

Let us estimate spinning down of magnetars in real, di- 
mensional units. We take as a base our main dimensionlcss 
parameters: angular velocity f2/Ox* = 0.7, Mach number 
M. = 2>, gravimagnetic parameter (3 = 10~^ and magnetic 
diffusivity fim = 10~^. 

Let us consider a neutron star with mass M = 1.4 Mq = 
2.8 X 10^^ g and radius Rns = lO'' cm. The density of the 
ambient interstellar matter is taken to be p = 1.7 x 10"^* 
g/cm^ (n = 1/ cm^). The sound speed in the ISM, Cs = 
30 km/s and a star moves in the interstellar medium with 
velocity u = 3cs = 90 km/s. 

Using the definitions /3 = 8ttp/Bq, p = pcl/j, we obtain 
the magnetic field at the surface of the "numerical" star 
Bo = {^TthPf'^p^'^Cs « 0.15 mcao G. 

We use the Bondi radius as a length scale. This radius 

is 

Rb = GM^/d « 2.1 X 10"cio^ cm . 

Recall that the size of the numerical star is _R, = 
0.025i?B ^ 5.2 X 10"c3,fcm in all simulations. Thus the 
magnetic moment of the rmmerical star is p = Borl/2 « 
1 X lO^^mc^Q Gcm^. Note that this is of the order of the 
magnetic moment of a real neutron star with surface field 
Bns ~ G and radius Rns = lO" cm. So we can suggest 
that a real neutron star is hidden inside numerical star. The 
neutron star with so strong magnetic field is considered as 
a magnetar. 
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The numerical star in our main case rotates with angu- 
lar velocity 0* = O.Tfix* = 2.5 x 10~* and this corre- 
sponds to a period of P« = 2n/Qr « 2.5 x 10^ s. 

For these parameters the corotation radius, Rcor = 
(GM/n.^)^/'^ = R,/uj,^^^ ^ 1.3R,, is appreciably less 
than Rm « 57?,. The light cylinder radius Rl = cP/2n = 
1.17x IO^^Pb cm is greater than the radius of magnetosphere 
taken from simulations Rm ^ 2.5 x lO^^c^Q^cm. So this case 
Rcor ^ Rm ^ Rl corresponds to the propeller regime. 

Therefore our simulations model an isolated magnetar 
which spins down by the propeller effect. 

Now we can estimate a time of spin down due to the 
propeller torque using equation (9). An angular momentum 
of a star is equal to L, = I,iU, where /, — 2/5AI,Ri - mo- 
ment of inertia for sphere. Taking /, « 10''^ gcm^ and il* = 
27r/P*, we obtain L, = 2-kI,/P, ^ 3 x 10"° Pg"^ gcm^/s. 

The angular momentum loss rate L = LLo, where for 
typical simulation run we find L ss 10"'^. Taking into ac- 
count dimensional parameter Lo « — 3x IO^^tiicJq' g(cm/s)^ 
we obtain L « —3 x lO^^nicj^Q* g(cm/s)^. 

In reality L is not constant and depends on many pa- 
rameters - velocity of the star, angular velocity and mag- 
netic field of the neutron star. We performed a number of 
simulations with difi^erent angular velocities ft* and differ- 
ent magnetic moments /x of the star and find dependencies: 
L oc M-^ " , or I, DC Pr' V°-®-M"°-*- Thus, 

L « 3 x 10^^nicJo'S?B«PB"'-^XJ°-* g(cm/s)^ . (17o) 

Thus the characteristic spin-down time is 

T = U/\L\ « 10* n^^d^B^t^Pl'Ml-^ yr . (176) 

Now we can estimate the time-scale of evolution at the pro- 
peller stage. For periods P* = 10^ s, which correspond to 
beginning of the propeller stage, the evolution scale will be 
AT = 10^ years, while at period P, = 10^ s corresponding 
to the end of propeller stage T = 3 x 10* years. Thus we 
sec that rnagnctars are expected to spin down very fast at 
the propeller stage. This time-scale however may be much 
larger if diffusivity is very small. 

5.1 Restrictions of the Model 

The axisymmetric resistive MHD numerical code used in 
this paper is a robust and has been checked at a number 
of tests and also it has been used for solution of similar, 
but somewhat simpler problems, such as Bondi accretion to 
a non-moving star, wind accretion to a moving magnetized 
non-rotating star, and others. This paper is a logical contin- 
uation of the series of simulations with increased complexity. 
The bow shock which forms in our simulations is similar to 
one obtained in simulations of the Earth's magnetosphere in- 
teracting with the solar wind (e.g., Gombosi, Powel, & van 
Leer 2000) and simulations of the pulsar wind interacting 
with magnetosphere of the companion neutron star (Arons 
et al. 2005). 

However, there is one factor which is important and 
basically is not known, this is the diffusivity. The role of 
diffusivity was underestimated in the analytical papers on 
propellers, however it is very important for mixing of the 
incoming matter to the magnetic field of the magnetosphere. 
In our model we took some value of diffusivity such that 



part of the incoming matter diffuses through the field lines 
and takes part in the spinning-down of the star. However, 
if the diffusivity will be much smaller, then the time-scale 
presented by the formulae (17b) will increase. 



6 DISCUSSION OF POSSIBLE MAGNETAR 
CANDIDATES 

It is not clear whether we can observe old, slowly rotat- 
ing magnetars. One of unknown aspects is that we do not 
know whether magnetic energy will continue to release, and 
at which rate, and what is the expected spectrum of such 
radiation. Theory of radiation from relatively fast rotating 
magnetars has been recently developed (e.g., Beloborodov & 
Thompson 2006), which predicts radiation from the whole 
magnetosphere mostly in the X-ray. However we do not know 
what would be the luminosity and spectrum of much slower 
rotating magnetars with periods P < 10^ seconds and thus 
it is not clear how to find them, even if many of them sur- 
vive. This is an obstacle in searching the old magnetars even 
with the best modern instruments. There are, however, can- 
didates for magnetars: SGRs and AXPs and also isolated 
neutron star (INS) candidate RX J1856.5-3754. 



6.1 Clustering of Periods of SGRs and AXPs 

One of amazing properties of SGRs and AXPs, which are 
strong candidates for magnetars, is the clustering of peri- 
ods around 6-12 seconds. This is still a stage when the light 
cylinder is smaller than magnetospheric radius, and the pro- 
cesses at the light cylinder may be responsible for the spin- 
down. The estimated age of these objects is around 10^ — 10* 
years which is an evidence of their fast spinning down at this 
stage. Analysis of possible evolution of periods have shown 
that after P > 10 seconds, these objects should either spin 
down much faster, or the magnetic field should decay very 
fast (Psaltis & Miller 2002, see also Colpi et al. 2000). The 
spinning-down at this stage has not been investigated yet. 
Recent relativistic 3D numerical simulations of aligned and 
misaligned rotators and further modeling in this direction 
may help to understand what is the expected spinning-down 
rate at the pulsar- type stages (Spitkovsky 2004; 2006). Un- 
fortunately, we can not extrapolate formulae (17b) to rrmcli 
shorter periods, P ~ 10 seconds, because our results arc only 
relevant to propeller stage with P > 10^ seconds. Our sim- 
ulations show, that propeller mechanism is very efficient in 
spinning-down of magnetars. However, we should remember, 
that there is a big uncertainty in the results of simulations: 
we do not know, which part of the ISM matter from the 
cross-section a = Trr^j will participate in the spinning-down 
of the star. This greatly depends on the diffusivity, which 
is not known. In spite of that it is important to know the 
dependencies of the torque on different parameters. So, our 
model can not explain the clustering, however if a magnetars 
with period P > 10^ seconds will be discovered in the future, 
then these formulae will describe the subsequent evolution. 



Spinning- Down of Moving Magnetars in the Propeller Regime 9 



6.2 Isolated Neutron Star candidate RX 
J1856. 5-3754. 

Of the small number of isolated neutron star candidates, 
the object RX J1856. 5-3754 is a special because it has a 
Ha nebulae with the shape of a shock wave (van Kerkwijk 
& Kulkarni 2001). The origin of this nebulae as well as the 
origin of the INS candidate are not yet known. One pos- 
sibihty is that it is a misaligned pulsar with magnetic field 
B ~ 10^^ - 10^^ G and period P Is. Then the nebulae may 
be a shock wave which appears as a result of the interaction 
of the relativistic wind of the pulsar with the ISM (e.g., van 
Kerkwijk and Kulkarni 2001; Romanova et al. 2001; Braje 
& Romani 2002; Romanova, Chulsky, & Lovelace 2005). A 
weak part of this model is the lack of pulsations. The X-ray 
pulsed fraction is < 1.3% (Burwitz et al. 2003). Another 
idea is that this is an old neutron star which had a recent 
accretion event and is slowly cooling down. Then the Ha 
nebulae may be connected with the ionization by the ISM 
(van Kerwijk & Kulkarni 2001). 

Another possible idea is that RX J1856. 5-3754 is a rela- 
tively old, slowly rotating magnetar (Romanova et al. 2001) 
which spun-down rapidly during its the propeller stage (Mori 
and Ruderman 2003). This model may explain the lack 
of pulsations. However, there is no simple explanation of 
the Ha nebulae. In spite of the large cross-section of the 
bow-shock, the energy of the ISM released at the shock 
Eshock ~ O-Spy'^TTrl, « 5 x lO^'^n^^^'vl^^ ergs s'^ is 
much smaller than that of the Ha nebulae. Also, reconnec- 
tion in the magnetotail (Toropina et al. 2001) gives insuf- 
ficient power. An additional heating mechanism of the star 
or magnetosphere is necessary to explain the Ha nebulae. 
Heating of the star due to accretion from the ISM is excluded 
because very little matter accretes to a strongly magnetized 
star owing to the very narrow accretion columns (Toropina 
et al. 2003). There is the possibility of heating of the star 
by the release of the magnetic energy (e.g., Beloborodov & 
Thompson 2006). However, the magnetic heating of slowly- 
rotating magnetars has not been investigated. Thus the pow- 
ering the Ha nebulae is uncertain but not ruled out. 



7 CONCLUSIONS 

Using axisymmetric MHD simulations we have studied the 
supersonic propagation through the ISM of magnetars in 
the propeller stage. We have done many simulation runs 
for the purpose of determining the angular momentum loss 
rate of the star due to the interaction of its magnetosphere 
with the shocked ISM. We conclude, that the interaction 
may be highly effective in spinning-down magnetars. A star 
with magnetic field B ~ 10^^ — 10^^ G is expected to spin- 
down in AT ~ 10'' — lO^years. This time may be longer if 
the ISM material does not efficiently interact with the ex- 
ternal regions of the magnetar's magnetosphere. Therefore, 
after relatively short stages of pulsar and propeller activity, 
a magnetar becomes a very slowly rotating object, with a pe- 
riod P > 10^ — lO^'s, which is much longer than the periods 
expected for ordinary pulsars. This may be a reason why the 
number of soft gamma repeaters, which are candidate mag- 
netars, is so small. We should note however, that the rate 
of spinning-down depends on the magnetic diffusivity which 



is not known. At lower diffusivity the rate of spinning-down 
will be lower. The INS candidate RX J1856. 5-3754 may be 
an example of a slowly rotating magnetar. However, this 
model does not explain the Ha nebulae. An ordinary mis- 
aligned pulsar explains the different features more easily, ex- 
cluding the fact that no periodic fiuctuations were observed 
from this object. 
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